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X(3872) and its production at hadron colliders 
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We evaluate the production cross sections of X(3872) at the Tevatron and LHC at NLO in a s 
within the framework of NRQCD factorization by assuming that the short distance production pro- 
ceeds dominantly through its x'ci component. The outcomes of the fits to the CMS pr distribution 
can well account for the CDF data, and are consistent with the value of k = Z C c-Br(X — > J /ip-K + -K~ ) 
constrained by the B meson decay data. In comparison, the molecule dominant production mecha- 
nism for X(3872) is inconsistent with both the pr distributions and the total cross sections of the 
CMS and CDF data. 



PACS numbers: 12.38.Bx, 13.25.Gv, 14.40.Pq 

Introduction- The hidden-charm state X(3872) was 
first discovered by Belle Collaboration in the J/ipTr + ir~ 
invariant mass spectrum of B + — > Jjjpir + it~ K + de- 



cay §, and confirmed by CDF [§, DO § and BaBar g 
Collaborations soon after. The closeness of the mass 
m x = 3871.68 ± 0.17 MeV § to the D°D*° threshold 
led many authors to speculate that the X is a D°D*° 
molecule g. (For recent reviews of X(3872), sec Rcf. 0) 

However, with the tiny binding energy Eb = {m D a + 
m B , a ) - m x = 0.142 ± 0.220 MeV §f, it is difficult 
to imagine that such a loosely bound state can have a 
large prompt production rate (i.e. not from B decays) 
in pp collision. In addition, DO Collaboration m found 
that the behavior of X (3872) production in pp collision 
is very similar to that of i/j', such as the cross section 
Pt distribution. In view of these, we proposed that the 
X(3872) is a mixed state of Xd(2P) (x' cl ) and D°D*° 
continuum pi . The two components both are substantial 
and may play different roles in the dynamics of X (3872): 
the short distance (the b- and hadro-) production and 
the quark annihilation decays of X proceed dominantly 
through the x' c i component; while the D°D*° component 
is mainly in charge of the hadronic decays of X(3872) into 
DDir, DD^j as well as J/ijjp and J/ipw. Based on the 

calculation for the ratio D , D ~ >Xcl „[ , we estimated that 

Br(B -* x'ci K ) = (2-4) x 10~ 4 g], which is consistent 
with the experimental constrains on Br(B — \ X(3872)K) 
renormalizcd by Z c5 , where Z c£ is the probability of the 
x'd component in the X(3872). We also calculated the 
X decays into J/ipp and J/tpui through rescattering of 
the intermediate states, i.e. the D°D*° component, and 
this may account for the observed large isospin viola- 
tion |Xo| . The closeness of the mass of X(3872) to the 
Z?°l)* fr threshold may be explained by the S-wave cou- 
pling between the J PC = 1 ++ cc component (x^i) an d 
the D°D*° component, which induces a sharp spectral 
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density and lowers the "bare" mass of x'd towards the 
D°D*° threshold |]J. 

As for the prompt production of X in pp collision, 
Suzuki [|12| also pointed out that the molecule assump- 
tion can hardly be consistent with the CDF measure- 
ment, based on an estimation of the wave function at the 
origin for the molecule. More explicit calculations given 
by Bignamini et al. |13| set the upper bound of the cross 
section of X as a D°D*° molecule to be 0.085 nb, which 
is smaller than the experimental lower bound 

a(pp -> X) x Br(X ->■ J/^tt+tt - ) = 3.1 ± 0.7 nb. (1) 

by about two orders of magnitude. 

Artoisenet and Braatcn K3 suggested that the DD* 
rescattering effects could enhance the molecule produc- 
tion cross section to values consistent with CDF data if 
one choose the upper bound of the relative momentum 
of DD* rescattering to be 3m v or more. With NRQCD 
factorization |D| at leading order (LO) and with the ma- 
trix element extracted from the CDF data of X(3872), 
they also made predictions for the molecule production 
cross section at LHC, which is about 3 to 4 times larger 
than the recent LHCb |6) and CMS JJ data. 

In our view, with the x'ci ~ D°D mixing model, 
one can naturally understand the prompt production of 
X(3872) at the Tevatron and LHC: the X(3872) produc- 
tion actually proceeds via the x'ci component and shares 
similar behavior to that of ip' . In this letter, we will study 
the prompt production of X(3872) at ncxt-to-lcading or- 
der (NLO) in NRQCD, and compare our results with the 
molecule model on both the total cross sections and pr 
distributions. We will also discuss the constrains on the 
value of Z C c indicated by our calculation. 

Prompt X(3872) production in NRQCD.-ln the x'd 
dominant production mechanism, the inclusive cross sec- 
tion of X in the J/?pir + ir~ mode can be expressed as 

da(pp -> X(J/4m + 7T-)) = da{pp -> X 'd) ' *, (2) 

where p is either a proton or an antiproton and k = 
Z c5 ■ Br with Br = Br(X -> J/tpTr+TT'). The feed- 
down contributions from higher charmonia (e.g. ip(3S) 



are negligible for the prompt production of V(3872)/Xci, 
so here "prompt" is almost equal to "direct", and the 
cross section of \' cl in (g) can be evaluated in NRQCD 
factorization, which is given by 



da(pp -> x'd) = E„ da{{cc) n ) 



(o„ 



(3) 



Sij.n / dxidx 2 Gi/pGj/pd&iij -> (cc) n )(0 



Xcl\ 



where Gij/ p are the parton distributions functions 
(PDFs) of p, and the indices i,j run over all the par- 
tonic species. The matrix element (On cl ) is marked by 
"n", which denotes the color, spin and angular momen- 
tum of the intermediate cc pair. Here we will evaluate 
the cross section at NLO in a s and LO in v (the relative 
velocity of cc in the rest frame of Xci)j therefore, only 



n = 3 P{ and 3 S{ are present here. 

Since Xci(l-P) and its radial excitation x'd share the 
same quantum number, they are physically alike. For 
the purpose of comparison, we follow the conventions of 
Ref. |18|, where similar calculations are done for Xci, to 
define the color-singlet matrix element as 



9 



(0^( 3 P 1 [1] )) = -|^ P (O)| 



47T 



(4) 



where R' 2 p(fi) is the derivative of the xj,i radial wave 
function at the origin, and the definition in (0) is different 
from that in |l5f by a factor of l/(2iV c ). Similarly, we 
parameterize the color-octet matrix element by the ratio 



r = m 2 c (O x '^ ( 3 sf ] ))/(O x '^ ( 3 PJ 1] )) 



(5) 



as in Rcf. |Tq|. Thus the cross section in (0) will be a 
simple function of parameters r, k and |-R 2P (0)| 2 . F° r 
further simplification, we fix the value of wavefunctions 



I^2p(0)| 2 Hi?U0)| 2 = 0.075 GeV 5 , 



(6) 



where in the second equality the value for Xc(lP) is cho- 
sen from the B-T type potential model calculation 119(1 
and was successfully used to evaluate the cross sections 
°f Xci.2 to account for the LHCb measurements 20 , and 
the first equality |i?2p(0)| 2 = l-^ip(0)| 2 i s assumed based 
on various potential model calculations |19[ , and this as- 
sumption has been adopted in our previous evaluation 

of the ratio D , D ^ Xcl ,a [PI . So it will be convenient to 

Br(B->-x<:iK) HJ 

compare the results of X production in pp collision with 
that in B decay. Actually, any change of |-R 2P (0)| 2 can 
be simply compensated by the corresponding change of 
k, and we will back to this point latter. 

As for the numerical calculation, we choose the same 
input parameters as Ref.|Q. we use the CTEQ6L1 and 
CTEQ6M PDFs |§l) for LO and NLO calculations re- 
spectively. The charm quark mass is set to be m c = 1.5 
GeV, meanwhile the renormalization, factorization and 
NRQCD scales arc /i r = p,f = \Jp T + 4m;! and [i\ = m c . 
To estimate theoretical uncertainties, we vary fi r and [if 



from mp/2 to 2top and choose m c = 1.5 ± 0.1 GeV. We 
refer other details of the calculations to Ref. jjj| . 

One should note that the D°D*° molecule can also 
be produced in pp collision through the short distance cc 
pair production with small relative velocity v. Therefore, 
the cross section can be factorized similar to that in (|3|) 
and given by |14| 

a{D°D*°) = &( 3 S [ i ] )(O dB ' ( 3 S l ^ ] ))+a( 3 sf ] )(0 DD ' ( 3 sf ] )), 

(7) 
where the PDFs have been integrated out in the co- 
efficients a( 3 S[ ), and the NRQCD velocity scaling 
rule p5[ and heavy quark spin symmetry p2[ have been 

adopted to truncate the formula to 3 S[ ' terms. The 
coefficients &( S S[ ' ) are the same as those in the char- 
monium cross section. We evaluate them at NLO in the 
same scheme mentioned above, and find that the ratio 
^SJ 11 )/^! 81 ) is about 5.3 x 10~ 4 for the CDF with 
p T > 5 GeV and 1.5 x 10~ 4 for the CMS with p T > 10 
GeV (the coefficient cr( 3 S[ ) has been divided by 2N C to 
match the convention in Rcf. WM). Note that the ma- 
trix elements in (R) are of the same order, and thus the 
color-singlet contributions can be neglected. 

Furthermore, the matrix elements (O dd ( 3 S[ ' ')) for 
a loosely bound state should be much smaller than those 
of charmonia, which can be justified by the calcula- 
tions of Refs. |13, 03], thus we neglect the contributions 
from (R) in our mixing model. However, Artoisenet and 
Braaten Q argued that the DD* / DD* rescattering ef- 
fects could enhance the molecule matrix elements to be 
consistent with the CDF data of total cross section for 
V(3872). Anyway, the two models are different in differ- 
ent combinations of channels in (||) and in (Q) . Thanks 
to the cross section pp distribution of V(3872) measured 
by CMS Collaboration [17]], which allows us to compare 
the two different combinations in the two models. 

Using formula (||) and (pi), we fit the CMS px data 
(VS = 7 TeV, \y\ < 1.2) §j\ by minimizing the % 2 , and 
the results are shown in Fig. 1 with the outcomes 



0.26 ±0.07, 



0.014 ±0.006, 



(8) 



where the central values correspond to % 2 /2 = 0.26, 
and the lager error-bars, indicated by the broad band in 
Fig. 1, are due to the insensitivity of the pt distribution 
to the parameter r in the range 0.20-0.40. Nevertheless, 
one can see that the central value of r in (0) is almost the 
same as that in Rcf. |||, i.e. r = 0.27, for Xd(l-P)- This 
strongly suggests that V(3872) be produced through its 
x'd component at short distances. In comparison, the 

Pt behavior of the 3 S[ ' channel is also shown solely in 
Fig. 1, which can hardly explain the data. Thus, the 
molecule production mechanism in (uft, where the con- 
tribution from the 3 S[ can be neglected as mentioned 
above, is disfavored by the CMS data. More explicitly, 
we use ([7|) to fit the pt distribution and get 



(O DD '( 3 S [ * ] ))Br 



6.0xlO~ 5 GcV 3 , x 2 /3=l-03. (9) 
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FIG. 1: The px distribution of prompt production cross 
section of X(3872). The CMS data are taken from Ref. 0. 
The theoretical curves are obtained from the two-parameter 
fits. The dashed, dot-dashed and solid lines represent the 
3 Sj , 3 P| and total contributions, respectively. The green 
bands denoted the uncertainties of the total results. 



For the CDF window (y/S = 1.96 TeV, \y\ < 0.6, p T > 
5 GeV), using the central values in (g), we predict the 
total cross section to be 



<r$ DF (pp^X(J/ilr. 



■)) = 2.5±0.7nb, 



(10) 



which is consistent with the data in (|l|). Besides, since 
the pt distributions of X(3872) and ip' production are 
very similar both for the CMS @ and DO data J§, one 
can expect that the same case would also occur for the 
CDF data. Therefore, we compare our prediction for the 
CDF p T distribution of X (3872) with that of V>' @, and 
the results are shown in Fig. 2, where the total cross sec- 
tion of ip' has been rescaled to be the central value in 
(pT|). As for the molecule production, not only the 3 S[ 
dominant mechanism can not account for the pt distri- 
bution of ip' as shown in Fig. 2, but also the predicted 
total cross section for CDF is too small, which can be 
obtained by using (R) and the matrix element in (M): 

(pp -r X( J/V>7r + 7r-)) = 1-1 ± 0-4 nb, (11) 



molecule 
a CDF 



which is about three times smaller than the CDF data 
in (0). This indicates that the CDF data and CMS data 
arc consistent with each other in our model, but seem to 
be inharmonious in the molecule model. 

For the LHCb window (VS = 7 TeV, 2.5 < y < 
4.5, 5 GeV < pt < 20 GeV), using the central values 
in (pi), we predict the total cross section to be 

VLHCbiPP "> X(J/^+tt-)) = 9.4 ± 2.2 nB, (12) 

which is about two times larger than the experimental 
data Il6|, 



°t%.(pP- 



X(J/^7r+Tr-)) = 5.4 ± 1.5 nB. 



(13) 



One should note that in (|T|) about 20% of the total cross 
section comes from b-decays, thus our prediction in (112]) 



FIG. 2: Comparison between our prediction and the CDF 
data tea] for the pr distribution. As mentioned in the text, 
the total cross section of the tfi' data has been rescaled to be 
the central value in (y_l|). 



is different from the LHCb data by about 2cr deviation. 
However, the error-bar in (O) is larger, and we expect 
more available data can be used to do the analysis. 

In fact, changing the values of r and k can improve our 
predictions, especially when the CMS pt distribution is 
not sensitive to r (as mentioned above). Thus, we fix r 
and fit k to the CMS pt distribution data, and the re- 
sults are shown in Table |, where only the central values 
of k and the predicted cross sections are listed. From 
the Table one can see that (Jcdf an d at LHCb nave dif- 
ferent r-dependence, and thus can not be simultaneously 
consistent with data in (Q) and (fl3|). 



TABLE I: The one parameter fit to the CMS pr distribution 
with fixed r. Only the central values of the obtained k and 
the predicted <y l CDF and u t hHCh are list here. 



r 


k 


xVd-o.f. 


ffcW(nB) 


crf HCb (nB) 


0.20 


0.021 


0.39 


3.26 


12.2 


0.25 


0.015 


0.17 


2.63 


9.87 


0.30 


0.012 


0.20 


2.28 


8.56 


0.35 


0.010 


0.27 


2.06 


7.72 


0.40 


0.008 


0.34 


1.90 


7.14 



The value of k can also be extracted from the B de- 
cay data, since in our model the branching ratio can be 
factorized as 

Br(B ->■ X(J/Tjm+Tr-)K) = Br(B -)■ x c \ K ) ' fc , ( 14 ) 

and the short distance branching ratio Br(B — > Xci-^0 
can be extracted by fitting the line shape of the experi- 
mental curves |24|, g£| . With a reasonable choice of the 
decay width of X(3872), the fit in Rcf. Q gives 

Br ftt (B -> x'eiK) = (3.7 - 5.7) x 10~ 4 , (15) 

which is consistent with our prediction in M. By com- 



paring (|lj) and (|l^) with 
Br{B -> X(J/i[>ir + Tr-)K) = (8.6 ±0.8) x 1(T 6 §, (16) 



the constrains on fc from B decay is given by 



A: = 0.018 ±0.004, 



(17) 



which is consistent with our fits to the CMS pr distri- 
bution for r = 0.20-0.26 in Table |. and thus consistent 
with the CDF data but not the LHCb data. 

As for the value of |i? 2 p(0)| 2 , if one choose a larger 
one than that in (|), say, 0.102 GcV 5 @, the value of 
k in Table Q will be decreased by a factor of 0.75 with 
fixed r, and the overlap between Table I and (Il7]) tends 
to disappear. Thus, our fits disfavor the lager one. 

Finally, with a modest value Btq = 0.05 which satisfies 
the experimental constraints [pi , the window of k in fll7| ) 
will correspond to the probability of x'ci component in 
the X(3872): 



Z c - C = (28-44)%, 



(18) 



which is consistent both with our original arguments and 
with the fits in Ref. §§. 

Summary. -Within the framework of NRQCD factor- 
ization, we evaluate the cross sections of A (3872) at the 
Tevatron and LHC at NLO in a s by assuming that the 



tion 

is almost 



short distance production proceeds dominantly through 
its x'ci component. The fit of the CMS pt distribu- 
data 1 17 gives the ratio r — 0.26 ± 0.07, which 
the same as that for Xci IP-Si ; an d strongly 
supports the x'ci dominant production mechanism for 
A (3872). The outcomes of the fits can well account for 
the CDF data [Q, [ll| , and are consistent with the value of 
k = Z C c ■ Btq constrained by the B decay data simultane- 
ously The predicted cross section for the LHCb is larger 
than the data |16| by a factor of 2, which might be due to 
the large theoretical and experimental uncertainties. We 
also evaluate the cross section of the D°D*° molecule |1J] 
at NLO, and find that the molecule dominant production 
mechanism for A (3972) is inconsistent with both the pt 
distributions and the total cross sections of the CMS and 
CDF data. 

We thank J.Z. Li, Y.Q. Ma and Y.J. Zhang for help- 
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Note : When this work was prepared, there appeared 
another paper p6j to study the prompt production of 
A (3872) as the Xci meson. Though some of their results 
are similar to ours, we put stress on the x'ci —D a D*° mix- 
ing model proposed in || , and find the x'ci component in 
A(3872) can be dominant if Br(X ->• J/^tt+tt") < 0.04. 
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